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Abstract 

We study the predictions of holographic QCD for various observable four-point 
quark flavour current-current correlators. The dual 5-dimensional bulk theory we con- 
sider is a SU (3) l x SU (3)r Yang-Mills theory in a slice of AdS§ spacetime with bound- 
aries. Particular UV and IR boundary conditions encode the spontaneous breaking of 
the dual 4D global chiral symmetry down to the SU(3)v subgroup. We explain in de- 
tail how to calculate the 4D four-point quark flavour current-current correlators using 
the 5D holographic theory, including interactions. We use these results to investigate 
predictions of holographic QCD for the AI = 1/2 rule for kaon decays and the Bk 
parameter. The results agree well in comparison with experimental data, with an accu- 
racy of 25% or better. The holographic theory automatically includes the contributions 
of the meson resonances to the four-point correlators. The correlators agree well in the 
low-momentum and high-momentum limit, in comparison with chiral perturbation the- 
ory and perturbative QCD results, respectively. 
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1 Introduction 



The success of the perturbative description of QCD allows us to understand the high en- 
ergy behaviour of strong interactions above 1.5 GeV. On the other hand, chiral perturbation 
theory (xPT) describes well the physics of strong interactions at low energy. In the inter- 
mediate region between both regimes, the situation is much less clear since neither of these 
theories behave perturbatively. One interesting and potentially powerful new idea to gain 
access to the non-perturbative regime of QCD is holographic QCD, which is based on the 
gauge/gravity duality [H El 13, H] . 

There are two kinds of holographic QCD dual models: there are 10D models based on string 
theory and supergravity [5] - [23], including studies of deep inelastic scattering [21] - [28J, 
and in addition, there are phenomenologically inspired 5D holographic dual models [29] - 
[33] . In both approaches, the description of confinement and chiral symmetry breaking has 
been tackled, and masses, decay constants, form factors and other properties of mesons have 
been calculated, yielding remarkably good agreement with experimental data. All of these 
estimates are based on two-point current correlators, which do not involve bulk interactions 
and pertain to the low-lying mesons. Given these initial successes, it is important that these 
holographic dual models of QCD are tested using processes that go beyond the properties 
of two-point current correlators, and include interactions in the bulk of the 5D theory. One 
such test is the computation in the 5D holographic theory of connected 4D four-point flavour 
current correlators, which can be compared with experiment and, in certain limits, with the 
results of chiral perturbation theory and perturbative QCD calculations. 

In this paper, we shall focus exclusively on such four-point flavour current correlators. 
These correlators are crucial to the resolution of a long-standing problem in QCD: the 
AI = 1/2 rule, which we will describe in detail in later sections and briefly here. In short, 
if one neglects CP-violating effects, there are two independent K° decays: K° — > 7r + 7r~ 
and K° —>■ tt tc . These two decays are combinations of AI = 1/2 and A J = 3/2 isospin 
amplitudes, Aq and A 2 respectively. Experimentally ReAo/Re^ = 22.2, and the largeness 
of this ratio is the AI = 1/2 rule. In the chiral limit, these two amplitudes are generally 
expressed in terms of the g 8 and g 2 7 parameters (see for example |36j), both of which depend 
on integrals over Euclidean momentum of certain four-current correlators. Our aim in this 
paper is to apply holographic QCD to calculate these observablesH 

The 4D theory we are trying to model is QCD with three massless quark flavours, possess- 
ing a global SU(3)l x SU(3)r symmetry which is spontaneously broken down to the vector 
subgroup via the quark condensate!! The AdS / CFT correspondence then immediately tells 

1 In a previous article [37], we briefly presented some of our initial results on this subject. In the present 
work, we significantly extend and improve upon our earlier study. 

2 Note that for simplicity throughout this paper we work in the chiral limit setting bare quark masses to 
zero, and ignore the anomalous and therefore explicitly broken C/(l) ax iai symmetry. We hope to return to 
these issues in a later publication. 
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us that the dual 5D theory should be a Yang-Mills theory with SU(3)l x SU(3)r gauge 
group, with a bi-fundamental bulk scalar field to provide breaking of this symmetry. The 
gauge fields in 5D couple to the QCD flavour currents, whereas the bulk scalar couples to 
the bilinear quark operator. In previous models, the inclusion of a bulk scalar field allowed 
a comprehensive description of chiral symmetry breaking [291 1301 l32| [33] . However, it is 
possible to take a limit of this theory where the entire description of chiral symmetry break- 
ing is encoded into the boundary conditions imposed on the gauge fields [HH [30] and the 
holographic theory contains only gauge fields in the bulk. This simplified holographic dual 
model turns out to be a reasonable approximation [3T1 ED], giving good results at least at 
the level of the two-point functions. The reason is that the condensate is an infrared (IR) 
effect, so that its influence can be modeled by an IR boundary condition. The complexity of 
the calculation of four-point current correlators in AdS / QCD means that this simpler form 
of holographic QCD, with only gauge fields in the bulk, provides an important starting point 
that can then be further refined. 

Our results are encouraging for AdS/QCD. As we discuss in detail in sections 5 and 6, 
we find that at leading order in a low-momentum expansion, the behaviour of the relevant 
correlators calculated in holographic QCD agrees with previous calculations using chiral 
perturbation theory, while at high momentum we obtain the behaviour predicted by per- 
turbative QCD. In the intermediate region, the momentum behaviour is governed by the 
exchange of meson resonances, and a significant advantage of the holographic calculation is 
that it automatically and consistently includes the contribution of the infinite tower of me- 
son resonances to the relevant correlators. Turning to a comparison with the experimental 
data, the results of a fit of the holographic predictions agree well, with an accuracy of 25% 
or better, which for quantities as difficult to calculate as the isospin amplitudes KeAo and 
Rev4 2 , is remarkable. Finally, we hope that the techniques developed here may be useful for 
more general calculations of n-point global symmetry current correlators in many AdS / CFT 
holographic dual models. 

This paper is organised as follows: In sections 2 and 3 we introduce the holographic QCD 
model that we use and present the relevant 5D propagators and interactions. In section 4 
we discuss the AI = 1/2 observables from the viewpoint of QCD, as well as a related but 
simpler observable Bk, parameterising K° — K° mixing. We shall also define the nature 
of the four-current correlators that we calculate using the holographic dual model, and the 
dependence of the parameters g$, #27 and Bk on these correlators. We also review the 
%PT predictions for the various observables. In section 5 we discuss the philosophy of the 
calculation and present the sum of the 5D Witten diagrams relevant for four-point functions. 
Section 6 contains our numerical results, and our conclusions are given in section 7, where 
we also briefly address the limitations of the model and possible avenues of improvement. 
Finally, four appendices contain technical details of the holographic calculation. 

Before we start upon our analysis we think it may be useful to offer the readers a "road 
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map" to follow the contents of this paper, depending on their particular interests. For readers 
interested in the AdS/QCD model and calculations of n-point correlators, and in particular 
four-point current correlators, sections 2, 3, and 5, supplemented with appendices A, B, C 
and D are recommended. For those interested in the physics of the kaon decays from xPT 
and perturbative QCD, section 4 is relevant. Readers interested in the comparison of the 
holographic calculation with the experimental results are directed towards section 6. 

2 The 5D Holographic Model 

Motivated by the AdS/CFT correspondence, and following on from the work of Refs. [291 1301 
I3"TI I3"2"l I3"3"] , we consider a 5D bulk theory defined in a constant curvature spacetime with 
the minimal field content as to describe current-current correlators in QCD. The spacetime 
metric is that of AdS^ space 

ds 2 = a 2 (z) (rf^ dx^dx u - dz 2 ) , (1) 

where a(z) = L/z, L being the curvature scale of the anti-de-Sitter space. The 5th- 
dimensional coordinate z holographically represents the energy scale of the 4D theory. We 
take z to extend from a UV boundary at z = Lq to an IR boundary at z = L\ > Lq. 

We are only interested in spin-1 4D operators such as qhl^qh and qRj^qR, where q 
can be u, d and s quarks. Using the well-known AdS/CFT relation between the dimension 
of such spin-1 boundary-theory operators A and the mass of the bulk vector fields m 5 , 
(A + 1)(A — 3) = m 2 , we find that A = 3 gives m§ = 0. We consider the chiral limit of 
QCD where the quarks are massless, so that global flavour currents are conserved and the 
boundary symmetry group is a global SU(3)l x SU(3)r. 

The rules of the holographic correspondence then tell us that the bulk theory is a pure 5D 
Yang-Mills with gauge group SU(3)l x SU(3)r. The boundary conditions on the UV brane 
z = L are such that the zero modes of the gauge fields in the \x directions are eliminated, 
so that no massless 4D gauge symmetry survives. The Lagrangian is given by [291 1301 13T] 

C 5 d = y/gM B Tr(-± L MN L MN - 1 R MN R MN ^j . (2) 

The scale M5 is some yet undetermined mass scale, g is the determinant of the metric and 
M = (/i, 5), where fi — 1, • • •, 4. The trace Tr is taken over the gauge group indices. 

We have Lm = L a M T a and similarly for Rm, where T a are the Hermitian generators for 
the Lie Algebra of the SU(3)l and SU(3)r groups, satisfying the following commutation 
relations and normalisations 

[T a ,T b ] =if abc T c and Tr[T a T b ] = 5 ab . (3) 
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The f abc 's are real and anti-symmetric in this basis. We write the following expressions for 
the gauge field strengths 

Lmn = 9mLn — 9nL m — i[L M , L N ] , (4) 
Rmn = 9mRn — 9nRm — i[Rm, Rn] , (5) 

which give us the following relation 

Lmn = 9mL n — OnL m + f abc L b M L c M , (6) 

and similarly for Rmn- 

We wish to work with the vector and axial-vector combinations of these gauge fields, so 
we define 

Vm = (Lm + Rm) , (7) 

Am = -^j= {Lm — Rm) ■ (8) 

The reason behind this choice is simple: the spontaneous breaking of chiral symmetry mixes 
the L M and R M gauge fields at the quadratic level. Therefore, the choice of basis as vector 
and axial-vector rather than left- and right-handed can be viewed as a diagonalisation of the 
equations of motion. Of course, at the cubic and quartic interaction level, there is mixing 
between Vm and Am, a fact which is integral to the calculation presented in this article. 

We can then express the Lagrangian above entirely in terms of vector and axial-vector 
fields. To eliminate the mixing between and V 5 and between A^ and A 5 , we need to 
include the following R^ gauge fixing terms 

%>f = -^Tr^d,V v -^d,{aV b )j , (9) 
C GF = -^Tr^d,A v -id,{aA^\ (10) 

where £ is the gauge parameter. 

We will go into unitary gauge in what follows, taking the limit £ — > oo at the appropriate 
stages of the calculations. Note that this will have different effects on the vector and axial- 
vector sectors, due to the different boundary conditions we impose, as explained below. In 
the next section, we present the full Lagrangian in terms of Vm and Am, and discuss the 
IR boundary conditions at z — Li, which correspond to spontaneously breaking the global 
SU(3)l x SU(3)r symmetry down to its SU(3)y subgroup. 
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3 Propagators and Interactions 



As we later discuss in detail in section 4, we wish to calculate certain four-point current 
correlators involving two left-handed and two right handed currents. We can expand the 
relevant desired correlator in terms of the vector and axial- vector currents Jy and J^, so 
that we can use the bulk Lagrangian in terms of the vector and axial-vector fields. 

According to the AdS/CFT correspondence, the boundary values of and are classical 
sources coupling to J v and J^, respectively. In order to calculate tree- level n-point functions 
for the currents, we need to solve the bulk equations of motion for the vector and axial- 
vector fields, substitute back into the action and treat this as the generating functional of 
the boundary theory. One thus has 



where, Sa<is is the Euclidean classical bulk action calculated with V^\uv = and A^uv = 
a M . We therefore need the bulk-to-bulk and bulk-to-boundary propagator for each of the 
gauge fields in the AdS field theory [21 EH] • The former allows us to construct the solution 
to the equations of motion from the interactions in the bulk of the AdS space, and the latter 
allows the construction of the solution of the equations of motion from the UV boundary 
value of the field. Green's second theorem gives us a straightforward relation between the 
two types of propagators. 

The procedure of finding the on-shell AdS action subject to certain boundary values of 
the fields can equivalently be formulated in terms of Witten diagrams, where one uses the 
vertices of the bulk theory to construct all the allowed Feynman diagrams connecting the 
boundary operators. The ingredients are the propagators and the vertices, as calculated 
from the bulk Lagrangian. In this section, we describe how to calculate the propagators. 
The vertices are simply derived from the full Lagrangian given in Eqs.f ll2[) - (TTTj) . We first 
justify our choices of boundary conditions for the various fields. 

Using the variational principle, the bulk equations of motion can be derived, along with 
the constraints that must be obeyed by any set of consistent boundary conditions. The 
UV boundary conditions on the bulk-to-bulk propagators can be chosen to be null Dirichlet 
for both the vector and axial-vector sectors. The IR boundary conditions distinguish the 
sectors, and allow chiral symmetry breaking (xSB) to be implemented into our model, by 
imposing null Neumann and null Dirichlet conditions on the vector and axial-vector sectors, 
respectively. This choice can be understood via an elegant argument: one can consider a 
bi-fundamental scalar living on the IR brane which acquires a vacuum expectation value. 
The effect of this on the boundary conditions is simple: it does not affect the vector sector, 
but changes the boundary conditions on the axial-vector fields from Neumann to mixed. 
This breaks the chiral symmetry spontaneously, and in theory we have one parameter to 
play with, analogous to the size of the quark condensate. Now, imagine removing the brane 



J d 4 xj!^(x)v l _ l (x)+J^(x)a l _ l (x)\ _ -S A dS 




(11) 
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scalar from the theory by allowing its mass to go to infinity The boundary condition on the 
axial-vector fields is now found to be a null Dirichlet condition. We also lose the parameter 
that allows us to tune the size of the symmetry breaking relative to the scale 1/L 1: which is 
set by the IR brane position. 

The final requirement is of course to account for the pions, which form a massless pseudo- 
scalar octet. We do this as follows: we impose warped Neumann boundary conditions on 
the A% field on both branes. This guarantees that even after going to unitary gauge in the 
axial sector, a zero mode remains in A 5 and cannot be gauged away. In the vector sector, 
we impose Dirichlet conditions on both branes for V5, so that going into unitary gauge here 
removes V5 from the theory. 



3.1 Interaction terms 

Here, we display the full boundary Lagrangian and bulk Lagrangian. Note that in the 
quadratic part, we have taken the limit £ — * 00 inside the differential operator acting on A u 
and Vft, but not in the A 5 operator. The quadratic part of the Lagrangian is used to calculate 
the propagators as shown below. The interaction Lagrangian provides the vertices needed 
to construct the Witten diagrams relevant for the calculation of any given n-point boundary 
current correlator. The full 5D Lagrangian can be written as a sum of the following terms, 
where the trace over the gauge indices is implicit 



■■boundary 



M 5 L 

2z 



1 Li 



vT (VvdsVp + A u d 5 A p - 2A p d u A 5 ) 



L 



(12) 



-■quadratic 



M 5 L 

2z 



(d 2 r)^ - zd z rT - Vff^ V u 



A p {d 2 ^ u - zd z tf - d^&^j A v 



+A 5 (-d 2 )A 5 + £A 5 d z [zd[-A 5 



(13) 



Uaa 5 ,va 5 a 5 = -*^§V" (d 5 V u [A 5 , A v \ + d 5 A p [A 5 , V u ] + d p A 5 [V v , A,]) , (14) 



'3— vector 1> 



M 5 L 
%/2z 



rf p r?° (d,V u [Ap, A,] + d,A v [V p , A a ] + d p A u [A P , K] + d p V u [V p , V a ]) , (15) 



--4— vector 



M 5 L 

Az 



rr-n va {V,y v [V p , V a ] + A P A U [A P , A a ] + AV p V u A p A a - 2V p V u A a A f _ 



+2V u A v V p A a - 2V p A v A a V p - A p V p A u V a - V p A p V u A a ), 



(16) 
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Cwmmaamm = "^V' (V„A s [V Vi A 5 ] + A,A,[A V , A,}) . (17) 

Examining the boundary Lagrangian in Eq. fTHZ]) . one immediately sees the familiar terms 
that are bilinear in and A^. These are responsible for the emission of vector and axial- 
vector resonances by the boundary theory current. The unusual term here is the mixing 
term between A^ and A 5 , which says that a boundary axial current can emit an A 5 particle. 
This term survives the application of the boundary conditions, and is in fact of paramount 
importance in the satisfaction of the Ward identities. This can be seen from a calculation of 
the axial current two-point function as described in appendix C. The Ward identity requires 
this correlator to be transverse, but this is only achieved by the AdS/CFT computation if 
one takes into account a diagram where the A§ field is emitted by one current and absorbed 
by the other. 

Finally, it is clear that the inclusion of the higher dimensional operators Tr(L\ IN + Rmn) 
and Tr(L A MN + R\ 1N ) in the Lagrangian of Eq.(j2J) results in other three-boson and four- 
boson vertices. However, one can show that they are sub-leading in the large M$L or, 
equivalently, the large N c limit, once one recognizes that M5L goes like N c parametrically, 
as in Refs. [251 I3U] . The argument goes as follows: the Lagrangian in five dimensions has the 
schematic form given by 



-5D 



'^Tr(F 2 ) + c ^Tr{F^ y/g~Tr(F 4 ) 



+ • 



where Tr(F 2 ) represents the leading terms of the gauge field theory as written in Eq. (E|) , i.e. 
Tt(LmnL mn + RmnR )i and the other terms signify higher dimensional operators, an 
example being Tr(g MS g N< ^g PR LMNLspLQFt + L — > R). Now, when we write this in terms 
of the fields A^, and A 5 , V5, we must remember that each factor of field strength F comes 
in with a factor of g MN . A factor of g MN brings with it a factor z 2 /L 2 . Recall also that 
^fg = (L/z) 5 . Thus, we can schematically write: 



1 z 

-Tr(F dd F dd ) + ci———Tr{F dd F dd F dd ) 

z (M 5 L) Z 



+C2jj-j~YTiTr{F dd F dd F dd F dd ) + 
{M 5 L) 



(19) 



where F dd simply means Lmn, Rain, i-e. the field strength with lower Lorentz indices. More 
generally, a gauge invariant operator with n factors of F will have a coefficient that goes 
like c n -2Z 2n ~ 5 / (M^L) 2( - n ~ 2 \ where all the c factors are of order one. This means that, in the 
large N c limit, the contribution from these operators is sub-leading to that from the term 
Tr(L MN L MN + R MN R MN ), as claimed. 
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3.2 SU(3)v sector propagators 

The vector bulk-to-bulk propagator has a transverse and longitudinal part, because we are 
working in the gauge V5 = (unitary gauge). The value of at the UV boundary is the 
classical source to which the 4D current Jy couples. On the IR boundary, we impose a 
Neumann condition on the vector field. We write (following |39|) 

(VV) = -iG v p {z, z') (rT ~ ^f) ~ z') {^pj . (20) 

Note that we are working in Lorentz indices, so we lower and raise these indices with the 
4D Poincare metric rf v . In Fourier space in the x M directions, but position space for the 5th 
direction, the propagators solve the equation 

dl ~ \d z + P 2 ) G^(z, z>) = J^6(z - z>) . (21) 

In addition, Gq (z, z') solves the same equation with p set to zero. The boundary conditions 
on Gp(z, z') are Dirichlet on the UV brane and Neumann on the IR brane, so that 

GjW)| 2=Lo = 0, (22) 
d z G^(z,z')\ z=Li =0. (23) 

From Green's second theorem, the bulk-to-boundary propagator is defined by the following 
limit 



(V»V V ) 

where 



dADS 



V) = _^L dz (y» V ») ; (24 ) 

Z Z=Ln 



The solutions are given by 



IT ZZ 

G^z, z>) z<z , = 2MsL{AD _ BC) i A MP^) + By^pz^CMpz 1 ) + Dy 1 (pz')) , (26) 



p V / I / 2 / U V I I 2 

dADS V V / V V 



z') = -iKj(z') IrT - ^ ) - iKUz') ) • (25) 



TtZZ 



G;(z,z') z>zI = ^YTTTlTFi ^[AJ 1 {pz')+By 1 {pz')][CJ 1 {pz)+Dy 1 {pz)\ , (27) 



where 



2M 5 L(AD - BC) 



A= -y^pLo), C = -y (pL 1 ), 

B= Ji(pL ), D = Jq(pL x ), (28) 



9 



and 

GftzyW = -^(z 2 - L*) , (29) 

G%{z,zr)^ = -^(J*-Ll). (30) 

Here J and 3^ are Bessel functions of the first and second kind in the conventions of Ref.[40j. 
From these bulk-to-bulk propagators, we find that the bulk-to-boundary propagators are 
given by 



W = ~L IAD-BC] ' (31) 
= 1. (32) 



Note that in calculating the bulk-to-boundary propagator, we use the bulk-to-bulk propaga- 
tor for z < z', so that 

K v {z >) = J^L dzG v^ z ') z<zl t (33) 

Z z=Lq 

(34) 



z=Lq 



3.3 SU(3)a sector propagators 

We here list the propagators for and A 5 . As explained above, the IR boundary conditions 
in this sector are chosen so that the SU (3) l X SU (3)r global chiral symmetry is spontaneously 
broken: 

Gf(z,z')\ z=Lo = 0, (35) 
G*(z,z')\ z=Li =0. (36) 

Similarly to the SU(3)v sector we define 

p^p u \ A , n / p^p u 



(A"A») = -iG*(z, z>) - ^) ~ iG£(z, z') ( ^- 1 . (37) 
The equation to be solved is the same as in the SU(3)y sector. The results are [29 | 130 1 I3T] 

Gf(z, z') z<z , = 2MsL( 7/- BC) [AJlipz) + fl Wl[^iW + Dy x {pz')\ , (38) 
z') z>zl = 2M ™^_ BG) WiW) + By 1 (pz')][CJ 1 (pz) + Dy 1 ( P z)) , (39) 
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where A and B are as for the SU(3)v sector, and 

C=-;Vi(pLi), D = J x {pL x ). (40) 

This gives us 

Gfc^-^-iS^I- (42) 



From these propagators, we obtain the bulk-to-boundary propagators as for the SU(3)v 
case, to find that 

K A ( r, Z' [CJ 1 { P z')+Dy i {pz')] 

K * {Z) = ~T [AD - BC] ' (43) 

tfW = Z -^4- (44) 



We also note that the current- current correlator as p — > is now given by (see appendix C) 

nx( P 2 )U = ^ = -^^. (45) 

This is a direct consequence of the IR boundary condition. The A 5 propagator in this gauge 
is simply given by 

(A s A s ) = -iG B p (z,*>), (46) 
where G^(z, z') is the limit as £ — > oo of the solution to the equation 

(^ - i-d z + + P ^ G 5 p (z, z') = -j^K* - A , (47) 
with the boundary conditions 

d z (aGl(z,z'))\ z=LQ = 0, (48) 

d z (aGl(z } z'))\ z=Li =0, (49) 

giving |f| 

G ^- M^m (?)• < 50 > 



3 Note the factor 2 difference between the A$ propagator here and the (incorrect) one used in our previous 
paper [57] . 
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4 Four-point Observables 



4.1 The AI = 1/2 rule 

Neglecting CP violation effects, there are two independent K — > tttt decay amplitudes, 
K° — > 7T + 7r _ and K° — > 7r°7r°. These amplitudes can be written in terms of the AI =1/2 
amplitude A and the AI = 3/2 amplitude A 2 as 

A(K ^ir + Tr-) = A e i5 ° + ^/T/2A 2 e i5 \ (51) 
A(K°^n°n°) = A e iS ° - V2A 2 e i5i . (52) 

The measured values of these amplitudes are 

ReA = 2.72 • 10" 4 MeV, ReA 2 = 1.22 • 1(T 5 MeV , (53) 

which gives 

I = 5^ = Re(jr -> (7T7r) J= o) = 

~ ReA 2 ~ Re(K -> (tttt)^) ' ' 1 J 

The large value of ReAo/Re^ is the so called AI = 1/2 rule. 

In the following, we use holographic QCD to calculate ReA and ReA 2 in the chiral limit. 
In this limit, at order p 2 in the chiral counting, all the AS* = 1 transitions can be obtained 
from the standard AS* = 1 effective Lagrangian, involving the usual g 8 and g 27 coupling 
constants (neglecting the small electromagnetic contribution, see for example [HJ 36J) 



y/2' 



£ eff =1 = —k V ^ V us [<fc£ 8 + 92 7 C 27 ] , (55) 



where 



with 



C 8 = and C 27 = ~(£ M ) 21 (£% + (C,) 23 (£") u , (56) 



i=l,2,3 



= -%Fl U{x)^D^U{x) , (57) 



and V u d = 0.974, V us = 0.224. The pion decay constant F n is taken in the chiral limit, 
where the masses of the u, d and s quarks are neglected (F n ~ 87 MeV). The matrix 
field U collects the Goldstone bosons of the spontaneously broken chiral symmetry of the 
QCD Lagrangian with three massless flavours, and D^U denotes the covariant derivative: 
D^U = d^U — ir^JJ + iUl^, in the presence of external chiral sources and r M of left- and 
right-handed currents. The parameters g$ and #27 encode the dynamics of the integrated- 
out degrees of freedom in the chiral limit. These include the heavy quark flavours as well 
as the light hadronic flavour states. Notice that the octet term proportional to g$ induces 
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pure AI = 1/2 transitions, while the term proportional to (727 induces both AI = 1/2 and 
AI = 3/2 transitions: 



Ao = -^=V ud V* s V2F n (g 8 + ^9wj ( M K ~ m l) > (58) 
A = -^V ud V: s 2F^g 27 (M 2 K -ml), (59) 

where Mr and are the masses of the kaon and the pion respectively, and Gp is the Fermi 
four-point interaction parameter. 

To calculate g 8 and #27 , we separate the long and short distance contributions as usual 
and perform an Operator Product Expansion (OPE), obtaining the effective Hamiltonian 
for \AS\ = 1 transitions [321 H3J Ej , 

G ^ 

^e// =1 = ~~r ^ u c *(A i )Q«(y u ) (/^ < m c = charm quark mass) , (60) 
* i=i 

a{fi) = Zi(fi) + Tyi(fi) , r = , £ q = V* s V qd . (61) 

The arbitrary renormalisation scale /i separates short- and long-distance contributions to 
the decay amplitudes. The Wilson coefficient functions contain all the information 

on heavy-mass scales. For CP conserving processes, the contribution involving the CKM 
elements of the top quark, encoded in yi(fi), is negligible and only the Zi(fi) are numerically 
relevant. The coefficient functions can be calculated for a scale \i > 1 GeV using perturbative 
renormalisation group techniques. They were computed in an extensive next-to-leading log- 
arithm analysis by two groups [45, 46J. After Fierz reordering, the local four-quark operators 
Qi(fi) can be written in terms of color singlet quark bilinears 

Qi = 4 s L ^d L u L ^^u L , Q 2 = 4 s L ^u L u L ^^d L , 

Q3 = 4 2j s L ^d L qLl^qL , Q4 = 4 ^ s L 7 M g L ?l7m^ > 

Q5 = 4 s L 7 M d L g R 7 M g H , Q & = -8 s L g fi , 

3 3 
<?7 = 4 ^ 2 e <? ^T^l 7ATf«SR ' Qs = -8^2 2 6q ^ L<lR ^ L ' ^ 
1 <? 

where the sum goes over the light flavors (q = u, d, s) and 

q R ,L = \{l ± l 5 )q , e, = (2/3, -1/3, -1/3) . (63) 
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The operators Qs, ■ ■ . ,Qe arise from QCD penguin diagrams involving a virtual W and a c 
or t quark, with gluons connecting the virtual heavy quark to light quarks. They transform 
as (8l,1r) under SU(3)l x SU(3)r and solely contribute to AI = 1/2 transitions. It is 
important to note that they are present only below the charm threshold, i.e. for \x < m c . 
Similarly the Wilson coefficients zj^ of the electroweak penguin operators Qj$ are non-zero 
only for /i < m c . Thus, in the following, only Qi and Q2 will be considered as we will always 
work in the regime \i > m c . Long-distance contributions to the amplitudes Aj are contained 
in the hadronic matrix elements of the four-quark operators, 

(Qi(p))i = (n7i,I\Q t (fi)\K°). (64) 

In the strict large N c limit, i.e. considering only the W exchange diagram with Z2 = 1 
we get g$ = #27 = 3/5, while experimentally, from Eq.(l53l and Eqs.( l5"8l -( l59l) one observes 
g^ xp = 5.1 and = 0.29. This shows how crucial the QCD dynamics is for the AI = 1/2 
rule. Important progress in the understanding of the AI = 1/2 rule was made when it 
was observed that the short- distance (quark) evolution, which is represented by the Wilson 
coefficient functions in the effective Hamiltonian of Eq. floTJl) . leads to both an enhancement 
of the 1 = and a suppression of the 1 = 2 final state. The octet enhancement [12] in 
the (QijQz) sector is dominated by the increase of z 2 when // evolves from Mw down to 
/i ~ 1 GeV, whereas the suppression of the A J = 3/2 transition results from a partial 
cancellation between the contributions from the Qi and Q2 operators. Taking into account 
the running of Z\ and z 2 between Mw and \i ~ m c = 1300 MeV, which gives z\ ~ —0.5, 
Z2 — 1.3, and still considering the matrix elements in the large N c limit, i.e. considering 
only the factorisable contribution, one gets #8 — 1 and #27 — 0.5. This gives values closer to 
the experimental ones but a factor 5 (3/5) is still missing for g 8 ((727)- Thus, perturbative 
QCD effects are far from sufficient to describe the AJ = 1/2 rule and QCD dynamics at low 
energies must be addressed beyond the leading N c limit, that is to say, at the level of the 
non- factorisable contribution [171 HH1 IIS EH [36] . 

Further progress was made when, in addition to the 0(p 2 ) weak AS = 1 Lagrangian 
of Eq. (1551) . the 0(p 4 ) AS = 1 Lagrangian was also considered. A full fit of all the weak 
Lagrangian constants was then carried out, taking into account not only the experimental 
K — > txtx amplitudes but also the experimental K — > 777777 amplitudes. It was found that 
the 0{p 2 ) contribution is expected to account for g 8 = 3.3 and g 2 7 = 0.23. The rest of the 
experimental amplitudes are expected to be explained by the 0(p A ) AS = 1 Lagrangian. 
Numerically, this 0(p 4 ) higher order (HO) contribution is equivalent to adding by hand a 
contribution of +1.8 and +0.06 to g$ and #27, respectively. In the following, we will only 
calculate the 0(p 2 ) Lagrangian constants g% and #27 > which should account for about two 
thirds of Re^o and three quarters of Re^. Therefore, for comparison with experiment, 
there are two equivalent possibilities: we either compare the values of g$ t 27 we obtain from 
Eqs.( l5"9"]) - (f70]) with the values 3.3 and 0.23 above, or, adding the 0(p 4 ) contribution by hand, 
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we compare the values we get for 

9l° T = 9s + flf °, (65) 
g$" = 9n + <&°, (66) 
with the values 5.1 and 0.29, with g>f° = 1.8 and g%° = 0.06. 

To calculate the non-factorisable contribution to #8 and (727, one can make use of the chiral 
symmetry properties of the AS = 1 effective Lagrangian of Eq. fl55l) . Instead of calculating 
the K — > mi amplitudes explicitly, it is much simpler to calculate them taking U — 1 in 
Eq.( l55i) . i.e. considering the processes with no external pseudoscalar and only two external 
sources coming from the covariant derivatives of U. For reasons explained in Ref.|36J it is 
convenient to consider the processes with two external right-handed sources (for instance 
for C 8 we consider C 8 3 J2i=i 2 3 -f^(^)2i (^ M )i3) • The non-factorisable contribution to this 
process of four-quark operators is then given by Green's functions involving, on the one 
hand, the two left-handed currents of the four-quark operator inducing this process, and, 
on the other hand, the two right-handed currents coupling to the right-handed sources. 
More precisely, including the leading N c non-factorisable contribution from Qi and Q2, the 
parameters g$ and #27 ar e given by the Q 2 integrals (with p the Minkowski momentum flowing 
between the two left-handed currents) of the two Green's functions [36, 15T]: 



WZS&iP) =£* * 3 / dVyd'z e^^ 2 )(0|T{L^(x)^( 2/ )L^(0X(z)}|0)| conn , 

(67) 

W^(p)=lim i 3 [ d'xdSd'z e ^^)(0|T{^ w (x)L^(0)i2| u (2/)^ s (^}iO}| conil) 



1^0 



(61 



with 



,2 



(69) 

'127 ( / ' ) = (■--]( / + Z 2 (fi) J ~gAS=2 (ji) , (70) 



and 



1 



A 4 



9AS=2(ji) = 1 - / dQ 2 W LRLR (Q 2 ), (71) 







while 



4Q 2 


f dQ p 


3 jp2 r l a l 3r lf iu / 


4/r 


1 Q 2 




3 pzHotfinia/ J 


47T 



LRLR\f)i 



(P), (72) 



) = -o^^^, -rrW^fCp) . (73) 
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Notice that we use the notation Q 2 for the Euclidean momentum (i.e. Q 2 = —p 2 ). 
From the above equations we see that the factorisable contributions to g 8 and g 27 are 

2 3 

9%{p) = -g^iO) + 5^2 (X), (74) 
3 

9$r(v) = g (^iO) + 2 2 (/i)) • (75) 

The remaining non-factorisable parts of Eqs. (|69[) - (|7ip are the subject of the calculations of 
this paper. 

In the above, we only consider, as necessary, the connected parts of the four-point func- 
tions. The currents are defined by Rq iq2 = gi7 M "V^^ Qz, Lq iq2 = ( ' 1 ~ 2 75 ' ) q2, and the sub- 
script of (?as=2 comes from the fact that this quantity also determines the AS = 2 transitions 
(see below) o 

In order to calculate the integrals of Eqs.(l69l)-(17TT). the Q 2 dependence of Wllrr and 
Wlrlr must be determined. However, this dependence is known only in the asymptotic 
regimes Q 2 — > and Q 2 — ► oo. In the limit when Q 2 — > 0, from x?T [52l I5T1 [36] . and after 
a long calculation, one gets 

W LRLR (Q 2 ) = 6 - 24(2/ x + 5/ 2 + Z 3 + l 9 )^ + (76) 

while for the other correlator only one group has calculated the ;\TT result [36J, 

W LLRR (Q 2 ) = -§ + (-y/3 + ^9) J| + (77) 

In these expressions the /j are the standard renormalized 0(p 4 ) chiral Lagrangian coefficients, 
usually denoted by Lj. 

In the limit Q 2 — > oo, and using Shifman-Vainshtein-Zakharov OPE techniques [43], one 
obtains [51, 36] 

o* F 2 

lim WW(Q 2 ) = +24vr 2 ^, (78) 

Q 2 ^oo 3 7T Q 2 3 7T iyQ 2 

where a s is the strong coupling constant, is the quark condensate, and l 5 is one of the 
0(p 4 ) chiral Lagrangian coefficients. Note that in Eq. (179|) the term depending on (ipip) is 
numerically dominant. Using 4D large- N c diagramatics, one can see that Wlrlr and Wllrr 



4 Notice that Eq. (|72[) has been modified by a normalisation factor 4, of which we were not aware in our 
previous paper [37] . This normalisation is correct when one sums over all possible planar flavour contractions, 
as we do here. 
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are given by a sum of simple to triple poles in Q 2 multiplied by polynomials in Q 2 . Combining 
this constraint with Eqs.(|76l)-(J771) and Eqs.(J78l ) -( J79l) . we then get the most general form for 
the Green's functions [SB EH] 



where the Mj's are the masses of the resonances and a i; 7$, and a-, f3[, 7^ are constants. 

Quite a few calculations have been proposed in the 1/N C expansion to estimate the Q 2 
integrals [I7J HEJ HH1 EH]- They all found a large enhancement of ReA together with a 
decrease of KeA 2 , so that the bulk of the AI = 1/2 rule can be explained. In those references, 
the size of the enhancement is determined essentially by two distinct factors. The first is 
the ;\TT behaviour at low scale, as determined by the chiral Lagrangian parameters, and 
the second is the size of the hadronic scales, namely the masses of the hadronic resonances, 
which will modify and terminate this behaviour at some higher scale. In particular, see 
Ref. [36J, which explains in detail why the interplay of the relevant hadronic scales — the 
small chiral constants and (tpip) on the factorisable side and larger resonance masses 
on the non-factorisable side — means that such a large non-factorisable contribution must 
be present. However, the effect of the resonances was not calculated explicitly in Refs. [T?] 
HH HHJ [36], but introduced in an indirect way. In [I7J 08] only xPT results (Eqs. (176 1 - ( 1771) ) 
were considered, with a cutoff put by hand at the mass of the resonances. The work of [36] 
used the form of EqB.flSOjl-flgll) to interpolate between Eqs. (j7Bjl - (j77)) and Eqs. (j75]l - (j75]l with 
a minimum number of resonances, while [49J employed Nambu-Jona-Lasinio models. The 
implicit assumption of these procedures is that the contribution to Wlrlr and Wllrr from 
the intermediate momentum region (0.5-2 GeV 2 ) can be obtained via a gentle interpolation 
(i.e. without "bumps") between the chiral behaviour and the OPE behaviour. In any case, it 
is clear that a method incorporating resonances explicitly would be highly preferable. This 
is precisely where the power of holographic QCD lies, at least for these observables, since 
it is a method where the effect of the entire tower of resonances for each channel can, in 
principle, be calculated. 

4.2 The Bk parameter 

From the Green's function Wlrlr-, there is another observable whose non-factorisable contri- 
bution can be calculated in the chiral limit and at leading N c order, and which can be used as 
a test of the holographic method described below. This is the observable, parameterising 
K° — K° mixing. At the quark level, K° and K° mix due to a box one-loop diagram where 




(80) 




(81) 
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the K transforms itself into a K through a pair of W bosons. This diagram leads to the 
following effective Hamiltonian |53j : 

C 2 M 2 

W3f =2 = K F ± + X " F2 + ZAcAtFs] C A5=2 (/i)g A5=2 (x) , (82) 

with 

<2as=2(» = (s L (x)7 M c/ L (x))(s L (x)7 Ai c/ L (x)), (83) 
and Cas=2 is the Wilson coefficient. From this effective Hamiltonian, defining 

< ^°|Qa5= 2 (0)|^ >= p 2 K M 2 K B K (fx), (84) 

the parameter Bk is defined as 

B K = C AS=2 (fi)B K (fi). (85) 

The large N c limit (i.e. the factorisable contribution) gives Bk = 3/4. 

In the chiral limit and at leading N c order, it turns out that the non-factorisable contri- 
bution is determined by the same integral of Wlrlr as the one found in parts of g$ and in 
5-27, Eqs. jESD-jZU). This gives 

3 

b k(h) = j9as=2{^)- (86) 

These relations come from a dynamical symmetry [M] relating part of the matrix elements 
of Qi and Q 2 with those of Qas=2- 

Unfortunately, there is no precise experimental determination of the Bk parameter. Thus, 
for our purposes, we will take Bk = 0.36 ±0.15 as a reference value, as obtained in the chiral 
limit in Ref . [36[ loTj [55] . Similar values have been obtained in the chiral limit, analytically 
in Refs.|56j and on the lattice in Refs.|57j. However, note that lattice calculations with 
physical quark masses [HE] have been shown to be sizeably larger than the chiral limit results, 
suggesting that the corrections beyond the chiral limit are large [221 EI] • 



5 Analytic Results 

5.1 Sum of the 5D Witten diagrams 

In this section, we show how to calculate the four-current correlators of Eqs.( l671) and ( |68|) 
in momentum space, i.e. 

WZAp) = i 3 lim(0|T{^(p)^| s (0^(-p)^(-0}|0) , (87) 

w£Trr(p) = ^ 3 |^(o|r{^ u (p)Z^(-p)Hi(0^(-0}io) • (88) 
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These expressions are in momentum space, so we use the tildes to refer to the Fourier 
transformed flavour currents. As we explain above, we use the vector and axial-vector field 
combinations, so that the vector and axial-vector currents are 

LUv) = ^(%» + r A , s M) , RUv) = ^(%» - J1,- S M) ■ (89) 

We have the following expansion for W^^ R {p) in terms of the vector and axial- vector currents 

£iim [j£(p) J v v {-p)~m jU-i) + j^p) Jv{-p)m jIh) - 
Jiip) Jv(-p)-m) jIh) - jkp) jy{-p)jm j?(-o - 

J v {p) J1{-p)JW) -JaH) - 4(p) Jl(-p)J2(i) JvH) + 

Jlip) Jl(-p)Jv(i) J?(-0 + Jlip) J u a(-p)J2(1) J p A {-l)\ , (90) 

and similarly for W^ U jf R (p). 

Having calculated in the previous sections the propagators for all the fields in our La- 
grangian, it is a lengthy but straightforward operation to construct all of the Witten diagrams 
for our four-point functions. One simply uses the bulk-to-boundary propagators to connect 
the four boundary points together through the vertices coming from the bulk interaction 
Lagrangian. Connecting points inside the bulk requires a bulk-to-bulk propagator. 

The inclusion of the boundary term Eq. fll2l) involving A 5 obviously increases the number 
of diagrams that contribute to any n-point function containing the axial- vector current. 
However, as shown in appendix D, the Ward identities satisfied by W^^rr an d W^xjj can 
be used to demonstrate that one gets the full result by considering purely the diagrams where 
only vectors and axial- vectors are connected to the boundary, and where only the transverse 
part of their bulk-to-boundary propagators is taken into account. This means that one need 
not consider the boundary term given by A^d^Ac, for the purpose of this paper. 

For both Wllrr(Q 2 ) and Wlrlr(Q 2 ), the 5D Witten diagram sum can be split into 
three distinct classes: diagrams where A 5 propagates in the bulk, X-diagrams involving the 
four-boson vertex, and Y-diagrams, which involve two three-boson vertices. Each class of 
diagrams contributes to the Green's functions at a different order of the momentum p: the 
A§ class contributes with order p° and higher, the X-diagrams to order p 2 and higher, and 
the Y-diagrams to order p 4 and higher. We refer the reader to appendix A for an example of 
each class of diagram, for the Jys d {p) Jy sd {l) J\ sdi~P) ^Asd( — contribution to the Wlrlr 
correlator]^ 



5 Note that we must respect the quark-flavour contractions, which eliminates some of the Witten diagrams. 
We then draw all the Witten diagrams which contribute to each term in this sum, and add all the various 
parts. It turns out that for Wllrr there are 36 distinct diagrams, which gets reduced to 24 diagrams upon 
enforcing the order of quark contraction. For Wlrlr, we find 40 such diagrams which give a non-vanishing 
contribution. These diagrams are the totality of planar diagrams when the order of quark flavour contractions 
is respected. 
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Once we sum the diagrams including all the contributions, we find that the two four-point 
functions are proportional to each other with a factor —16. This factor comes from the 
£77(3) group theory structure. The proportionality is strictly correct only in the l a — > 
limit, which is the limit required for the computation of the g% and #27 parameters. We 
therefore have 



W LRLR (Q 2 ) 



4iQ* 
3 Fl 



E(p = iQ), 



W, 



LLRR 



(Q 2 ) 



±91 

12 Fl 



E(p = iQ), 



(91) 
(92) 



where £ denotes the sum of the diagrams and can be written as E = Ex+E^+Ey, referring 
to the distinct classes of diagrams. The Ex and Ea 5 components are given by 



Ex(p) 
and 



M 5 L^ [ d - l] 3 /"dz 
d 



[< 2 + < 2 ] [< 2 + Kf] - 4<^«) , (93) 



V2 J d 

where d = 4 is the dimension of spacetime and 

A'(z,z') = 2(K^d z K^-K v d z K p A )(K^d z ,K^-K^d z ,K^ 
+ (KfaK* - K v p d z K*) {K*d z ,K v p - K v p d z ,K£) 
+ {KfaKZ - Kq8 z K p ) (KfaKX - K v d z ,K*) . 



(94) 



(95) 



As for the Y-diagrams, the integrations are more involved, but the sum can be written as 
E Y (p) = Y 1 + Y 2 + Y 3 + Y A - 2Y 5 - 2Y 6 , where 



Y 2 
Y 3 
Y 4 
Y 5 
Y e 



—1 



—1 



M 5 L \ 

V2J 
M 5 L 

V2 
M 5 L 

V2 
M 5 L 

V2 
M 5 L 

V2 
M 5 L 



[d-l]p 2 J y / ^(z)G^(z,z')K^z')K v (z)K^z 
[d~l]v 2 f y / ^KV{z)G*{z,z')KV{z')K${z)K${z 

[d~l]p 2 J y J -^-Kp(z)Go(z, z')K p a (z')Kq(z)Kq(z' 
[d-l}p 2 J y J ^K^z)G^z,z')K^z')K v (z)K^z 
[d~l]p 2 j y / ^K^z)G^(z,z')K^(z')K^z)K v (z 



(96) 
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Figure 1: A plot of Wlrlr against Q, in the limit L — > 0. The correlator Wllrr is —16 
times smaller. 



We perform all the integrals with the limits L\ and L , and show the full results in appendix 
B. The results are very complicated expressions, but one can then take the limit Lq — > 
smoothly. All the divergent contributions cancel, and we obtain the simple result 
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(97) 



where Jo,i = Io,i(QLi) are the modified Bessel functions of zeroth and first order, respectively, 
and Q is the Euclidean momentum. This simplified form is more appropriate for the analysis 
of the high and low Q behaviour of the correlators. See Figured] for a plot of Wlrlr against 
momentum, with 1/L± = 280 MeV. Note that Wlrlr is found to be positive definite, while 
Wllrr, is negative definite because of the proportionality. Both correlators also approach 
zero as Q — ► oo, and satisfy the "sum of poles" functional form of Eqs. (!80]) - (!5T]) . More 
precisely, the high Q behaviour of the correlators is given by 1/Q 2 , which is the correct 
functional form predicted by perturbative QCD, Eqs. (1781 - (1791 . 



5.2 The limit Q — > and connection with chiral perturbation the- 
ory 

The pole structure of the propagators for low momentum constitutes a strong check on our 
calculation. Another check is whether our results agree with x?T which, as explained above, 
gives us a constraint on the behaviour of the correlators as Q — ► 0, Eqs. (I7B1) - (1771) . Taking 



21 



that limit in the expression of we obtain 

Jto E(Q) = W (gL + f - |^L? + 0&)) . (98) 

This is indeed the functional form required by ;\RT, the Q 2 pole being due to the massless 
pions. Our correlators therefore have the low Q behaviour given by 



lim W LRLR (Q 2 ) = 6-^^ + 0(Q 4 ), (99) 



limHWQ 2 ) - -§ + ^g + (Q«). (loo) 

This is to be compared with the expressions obtained via xPT in the chiral limit, Eqs. flTS"] )- 
f!77|) . A plot of our results versus those of x?T makes things clearer, for a value of 1/Li = 280 
MeV (Figj2]). The matching obtained for Wlrlr is very good for the range of validity of 
%PT, while Wllrr does not exhibit as good a matching (see below). Note also that the xPT 
results shown in the plots do not contain any 0(p e ) contribution, while our 5D result is to 
full order in p. 

In Ref. [21], the coefficients of the 0(p A ) chiral Lagrangian were calculated in an AdS 
setting with identical field content to the one used in this work. This allows us to compare 
our predictions in the low momentum limit to those of %PT with the AdS U coefficients 
calculated in Ref. [3T] • Using those results, and the relations found between the /, coefficients, 
i.e. I2 = 2/i and l 3 = —6/1, Eqs.(l99 l -f ll00p can be rewritten as 

O 2 

Jim W LRLR {Q 2 ) = 6 - 24(2/x + U 2 + h + h)j^ + - (101) 



and 

lim W LLRR (Q 2 ) = - 3 - + h 9 - k)% + - (102) 

Notice that the first expression coincides exactly with the pure calculation, Eq.f lTBl . 
Similarly, for Wllrr, the 0{p 2 ) coefficient —3/8 and the 0(p 4 ) l 9 coefficient +3/2 coincide 
with the corresponding coefficients of Eq.( l77l) . However, the holographic calculation does 
not reproduce the 0(p 4 ) l 3 coefficient of Eq. (!771) . yielding a factor —3/2 in place of —15/2, 
so that the total Q 2 coefficient in Wllrr differs by a factor two approximately from the xPT 
result. We do not understand this discrepancy. We have performed a variety of consistency 
checks on the 5D calculation, and we do not see any possibility of deviations which would 
alter the proportionality between Wllrr and Wlrlr- This makes us confident that our 
results are correct. It seems possible to us that the problem might lie with the sole and 
rather subtle x?T calculation of the / 3 dependence of Wllrr- Note also that this difference 
in the Z 3 coefficient for Wllrr is not significant enough to alter the fact that we find a large 
enhancement for g$ below. 
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6 Numerical Results and Discussion 



In this section we present our numerical results for the AI = 1/2 rule and the Bk mixing 
parameter. To obtain the values of the parameters g$ and g 27 , we must integrate the two 
four-current correlators over the Euclidean momentum as explained in Eqs. (|69p - (|7ip . This 
integral should ideally be regularised in the same scheme as the Wilson coefficients of the 
four-quark operators responsible for the kaon decay. These are z\ and z 2 , and they are 
usually renormalised using dimensional regularisation in three distinct schemes: Leading 
Order (LO), t'Hooft-Veltman (HV) and Naive Dimensional Regularisation (NDR). On the 
other hand, the sharp UV boundary at z = Lq in the 5D calculation implies that a hard 
cut-off at scale 1/Lq should be employed in the momentum integral, and that one should 
use the appropriate values of the Wilson coefficients z\ and z 2 at this energy scale. We 
therefore focus upon the Wilson coefficients calculated in the LO renormalisation scheme as 
this provides a closer, though admittedly not exact, match to the holographic part of the 
calculation. One may justify this choice as follows: our results for the correlators have the 
functional form expected from QCD calculations in Eqs. ( 1781) -(1791). The integrals carried out 
in the computation of g$ and g 2 i will therefore have a logarithmic divergence with respect 
to the cut-off, rendering the integration stable under changes of the high-momentum scale. 

We employ two choices of the high-momentum cutoff: 1300 MeV, which is approximately 
the mass of the charm quark, and 1500 Mevj§ Below, we also demonstrate that the results of 
our calculation are stable against changes in this cutoff. The values of the Wilson coefficients 
z\ and z 2 that we use, as calculated in the work of [HO] for the LO scheme, are Z\ = —0.625, 
z 2 = 1.345, at 1300 MeV and z 1 = -0.5699, z 2 = 1.307, at 1500 MeV. For C AS =2 we use the 
values 1.17 and 1.21, respectively. 

Below, we employ a self-consistent prescription to carry out two distinct fits. In the first, 
we fit to the following observables: the mass of the rho meson m p , the mass of the a± axial- 
vector meson m ai , F n , g 8 and g 27 . The values of F n , m e p xp and m^ p that we fit to are 87 
MeV, 776 MeV and 1230 MeV, respectively. For g 8 and g 27 , as explained in section 4, the 
values we fit to are g 8 = 3.3 and g 27 = 0.23, equivalent to gJ OT — g% + g 8 ° = 5.1, and 
#27° T = 927 + 927° = 0.29, Eqs. (l6"5"]) - (I6"6"]) . In the second type of fit, we fit to u> and B K 
instead of g 8 and #27, taking the values to = 1/22.2 and Bk = 0.36, see section 4. The only 
free parameters in both fits are M 5 L and L\. Inside the integral, we use our full result for 
finite L Q , the Minkowski version of which is shown in appendix B. The predictions of the 

6 In taking the cutoff to be 1500 MeV, which is above the charm threshold, we should also consider the 
contribution of four-quark operators involving the charm quark (see. e.g. [60]). However, for a scale not 
larger than 1500 MeV, which is well below the mass of the charmed resonances, their contribution is expected 
to be quite small, and in the following we neglect them. 
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model for F n , m p and m ai , entering in the fits, are 



F *=^T^k' (103) 



and, to a good approximation in the range of interest, 



th 2.12(^-0.282^) 
m P « -j— r\ ' ( 104 ) 



m ^ (L 0.085L„, 



Observable 


A 


B 


C 


D 


m ai / m aT P 

„TOT / „ e:E P 
i/8 /i/8 
„TOT I n ex P 
i/27 / i/27 

1/w 

TDth 

K 


1300 MeV 
274 MeV 

0.91 

0.95 

1.15 

0.74 

0.85 

19.5 

0.38 


1500 MeV 
275 MeV 

0.90 

0.93 

1.17 

0.72 

0.85 

19.2 

0.38 


1300 MeV 
277 MeV 

0.93 

0.97 

1.12 

0.75 

0.79 

21.4 

0.34 


1500 MeV 
280 MeV 

0.92 

0.95 

1.14 

0.74 

0.78 

21.3 

0.34 



Table 1: Columns A,B show a fitting to m p , m ai , F^, g& and #27- Columns C,D show a 
fitting to m p , m ai , F n , 1/u and B K . Note that, as explained in section 4, gJ OT — g$ + 1.8, 
927° T = 927 + 0.06, where g$ and (727 are the quantities we calculated from the AdS model. 
We use the values F n = 87 MeV, mf p = 776 MeV, r<f = 1230 MeV, g e 8 xp = 5.1 and 
927 P = 0.29. 



The results for both fits are given in Table 1 and are quite similar. Taking into account the 
relative crudeness of our model, and the use of the large N c expansion of QCD, we find these 
results quite good. The discrepancy with experiment never exceeds ~ 25%, and for some 
observables is much smaller. It must be emphasised again that, having picked the values for 
the upper cut-off to be 1300 MeV and 1500 MeV and thereby fixed L Q , the only remaining 
free parameters are L\ and M 5 L. Therefore, we fit five independent observables with only 
two free parameters. 

To understand the structure of the g 8 and g 2 7 results in Table 1, it is useful to decompose 
the numerical results into three components: the leading N c factorised part of Eqs.fl69l ) -( l7"Ti) 
given explicitly in Eqs. (!74|) - (!75|) . the 1/N C non-factorisable contribution of Eqs. (!69!) - (!7T1) 
(encoded in the Wlrlb, an d Wllrr Green's functions we calculated), and the contribu- 
tion from the higher-order corrections, gf° and g^ , which as said above we take from 
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Ref. [41]. For example, for the results of column A this decomposition goes as follows: 
gJOT = L06 + 92 + i 8 = 3 78 and g TOT = QA3 _ o 24 + .06 = 0.25, whereas the 

experimental values to compare with are g^ xp = 5.1 and g^j p = 0.29, respectively. We see 
that, for g$, the non-factorisable chiral limit contribution of 0.92 is of the same order as 
the factorised contribution of 1.06, and effectively doubles it. For #27, the non-factorisable 
contribution of —0.24 effectively divides the factorised result of 0.43 by more than two. As 
already mentioned above, the fact that the 1/N C contribution can be as large as the factorised 
part, even though the 1/N C series is still expected to converge, can be naturally explained 
by the interplay of the various hadronic scales involved in g 8 and $27 |36j . 

For a rather conservative estimate of the erroi0 involved in our calculation, it is interesting 
to compare the results of Table 1 with the values we obtain by fitting only F w , m p and m ai . 
In this case, for example with L„ 1 = 1300 MeV, we get FfjF^ v = 1.00, m p h /m e p xp = 0.98, 
m^Jm^P = 1.02, and gJ OT = 1.06 + 1.32 + 1.8 = 4.18, which gives gJ OT /g e 8 xp = 0.82. In this 
case, the non-factorisable chiral limit contribution of 1.32, is ~ 40% larger than the value 
0.92 above. For g 2 7 we get a smaller result: g\° T = 0.43 — 0.35 + 0.06 = 0.14, so that the 
non-factorisable chiral limit contribution of —0.35, is also ~ 40% larger than the value —0.24 
above. The small total result it gives is not surprising, as #27 involves the difference of two 
large positive contributions, i.e. a 1/N 2 correction of order ~ 30% of our 1/N C contribution 
would bring g 2 -j close to the experimental value. 

Note that our results are similar to what has been obtained in other analytical calculations 
utilising the 1/N C expansion [HI [361 [49]. I n particular, in Ref. [36] the values g1° T /g^ = 0.76 
and 5 , 2 Z 7 OT /5 l 27 P = 0-^9 have been obtained. The two methods are, however, quite different, 
as explained in section 4. The main advantage of our model is that it allows the calculation 
of the four-point functions in the entire relevant momentum region within one consistent 
setting, thereby removing the need for interpolation in any specific momentum range. 

One could ask why we took the parameter M 5 L as a free parameter in the fits. In Refs. [611 
[2nHS01EIlE21E31EZ], it has been shown that this parameter determines the high Q 2 logarithm 
of the vector-current two-point correlator, so that it can be fixed from a matching with the 
corresponding QCD coefficient [62], yielding 

M 5 L = . (106) 

The point is that, for two-point functions only, there is enough parameter freedom to match 
the QCD logarithm (see in particular [2211301 EI]), unlike the more complicated case presented 
here. Clearly, it is not expected that the model we consider, based on a simple slice of AdS 
with a hard cut-off in the UV at the scale Lq 1 , would lead to the exactly correct QCD 
behaviour. However, comparing our results with those of QCD in Eqs. ( 178"]) - ([79"]) . we observe 
that for both Wlrlr and Wllrr we reproduce the good functional behaviour (i.e. the l/Q 2 

7 This is also useful as an estimate of higher order 1/N C correction effects. 
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dependence). We also observe that we get the correct sign for the coefficients of 1/Q 2 and, 
amazingly, we even get values for these coefficients which are within a factor 2 — 3 of the 
perturbative QCD result, for Q 2 ~ 2 — 5 GeV 2 . Note that this is a very surprising outcome, 
especially for Wllrr, because the high-momentum dependence of the latter involves the 
quark condensate contribution, as seen in Eq. ( 179ft . In fact, the dominance of the quark 
condensate term in Eq. (179"]) guarantees that Wllrr is negative in the far UV, and so we 
predict its sign correctly, although our model has no equivalent of the quark condensate (or 
of a s for that matter). As mentioned above, one way of introducing a tunable condensate is 
by adding a scalar field in the bulk. It would be interesting to see whether the calculation 
of Wllrr in that case gives a more accurate description of the high- momentum behaviour. 

We showed in section 13.11 that the 5D bulk Lagrangian we employ is the leading order 
Lagrangian in the large-iV c expansion, and that operators of higher mass dimension are 
sub-leading in N c . In principle, these operators may contribute to the four-point functions 
calculated here. In this section, we have presented the results of a fit of five observables using 
only two independent parameters, which are the IR brane position L\ and the dimensionless 
combination M^L. Thus, the fact that a fit using only the leading operator of Eq.(j2]) gives 
good agreement to the data is non-trivial. If we were to introduce the full set of sub-leading 
operators that contribute, then we increase the number of free parameters of the model 
(because the coefficients of the new operators are unconstrained by bulk gauge invariance or 
other symmetries) and the fit loses predictivity. Therefore, the success of the restricted fit 
performed here shows that the coefficients of the sub-leading operators are not anomalously 
large. 

A technical, but important, issue to mention concerns the gauge symmetry of the 5D 
theory. To carry out the 5D calculations, we had to choose a specific gauge to work in, and 
we picked the convenient gauge taken in the limit £ — > oo. Now, it is clear that the 
results of the holographic calculation must be independent of the gauge parameter £, and of 
any choice of gauge- fixing. This must, in fact, be a feature of any holographic calculation 
involving gauge freedom in the AdS theory. In previous AdS / QCD computations there had 
been no need for concern, since these calculations only involved two-point functions and all 
the propagators in the 5D theory were boundary-to-boundary ones. Boundary theory current 
conservation then automatically projects out the longitudinal gauge-dependent component of 
the 5D propagators, rendering the results gauge-independent. Unfortunately, the situation is 
less clear in the case of four-current correlators, because the latter involve Witten diagrams 
with explicitly gauge- dependent propagators (exchange diagrams |38j). In the above, we 
trust that the power of AdS/CFT guarantees that any such holographic calculation will 
yield results that are independent of the 5D gauge fixing. 

The model has many shortcomings, due to its crude nature. For example, a concern for 
this class of models is the behaviour of the masses of the resonances My n and M^ n as n 
approaches infinity. One finds that the simple treatment presented here shows that the 
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masses of the Kaluza-Klein modes go like n for large n, in sharp contrast to the predictions 
of large-iV c theories, in which the masses of the resonances go like \fn. A recent paper 
[U3] has shown, however, that with a more sophisticated handling of the IR truncation of 
the AdS space, one can indeed recover the large- N c Regge behaviour. Whether this will 
improve the results obtained here remains to be seen. Finally, one must keep in mind that 
the calculations done here were all in the chiral limit. In order to account for massive quarks, 
one would have to introduce the bi-fundamental scalar of Refs. [29], ED] ■ In particular, this 
would allow the inclusion of the mass of the strange quark, whose effect could easily be 
as large as 20-25% for g 8 , or even more for g 2 7, due to the cancellation of factorizable and 
non-factorizable contributions in the latter. 

7 Conclusions 

In this paper we have calculated, within the simplest possible version of holographic QCD, 
the four-point flavour current correlators crucial to the resolution of the AI = 1/2 puzzle of 
QCD. We believe that our results are quite encouraging for the AdS/QCD approach. The 
holographic theory automatically and consistently includes the contributions of the infinite 
tower of meson resonances to the four-point correlators. We also reproduce, to a good level 
of accuracy, the low-momentum and high-momentum behaviour of these correlators, as de- 
duced from chiral perturbation theory and perturbative QCD, respectively. This agreement 
is particularly impressive for the correlator Wlrlr- Moreover, the results of a fit of the 
holographic predictions to the experimental data agree well, with 25% accuracy or better, 
showing that the dynamics of the A = 1/2 rule is operative in AdS/QCD. For quantities 
as difficult to calculate as the isospin amplitudes of kaon decay ReA and ReA 2 , this is 
remarkable. 

A rather obvious limitation of the model concerns the description of x^B. As explained 
above, although the imposition of IR boundary conditions on the bulk SU(3)l x SU(3)ji 
gauge fields correctly incorporates the leading xSB behaviour, a bi-fundamental bulk scalar 
is needed to fully account for the physics of xSB. The inclusion of this field will directly 
introduce pseudo-scalar resonances into the 4D field content, and these will indeed have 
relevant contributions to the four-current correlators calculated here. We will also have 
an extra parameter that can be tuned [29], [30], corresponding to the quark condensate. We 
have also not included the effects of the anomalous U(1)a symmetry of QCD, nor the explicit 
breaking of chiral symmetry due to bare quark masses. One envisions these improvements 
having a complicated yet positive effect on the calculation of four-point current correlators 
presented in this paper. 

We believe that the results of this paper for four-point functions show that it is worth 
investigating further the predictions of AdS/QCD. 
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Appendices 



A The different classes of 5D Witten diagrams 

We adopt the following conventions in labeling the vector, axial-vector and pseudoscalar 
propagators: 

z z < — (A^A") 

Using the vertices calculated from the full Lagrangian, we can then draw all the Feynman 
diagrams that contribute to a specific four-point function in momentum space. To illustrate 
this, we consider here the < Jy.ly.1 aJ a > piece contributing to Wlrlr- This is given by the 
following sum of diagrams: 



sd, jj, ds,a sd,jj, ds,a 




ds,/3 sd,v ds,/3 sd,v 
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In principle, one could have two more "cross diagrams" , but they are eliminated by the large 
N c condition that diagrams must be planar. 
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B The sum of Witten diagrams for finite L 



o 



Here we present the results of the Witten diagram summation for finite Lq in Minkowski 
space. As in section 5, we write E = Sx + £U 8 + Sy. Defining 

C n {z) = CJ n (pz) + Dy n (pz), 

c n (z) = cj n (pz) + Dy n ( P z), 

F(z) = p 2 [z 2 C 2 (z) + L\C {z)l 

G{z) = p 4 z 4 C 2 (z)-2p 3 z 3 C 3 (z)-p 4 L 2 lZ 2 C 2 (z). 



we obtain 
and 

S x = 3iM 5 ( - 
1 



36iM 5 L 



Ms 



L 2 {L\ - LIY{AD - BC) V 



[z 2 C 2 {z)] 



Li 



L 



-{C^zf - C (z)C 2 (z)) 



+ 



+ 



L 2 (AD - BC) 2 \ 2 

^(z) 2 - C (z)C 2 (z)) 



Lq(AD - BC) 2 V 2 
1 

L 2 (L 2 -L 2 ) 2 (AD-BC) 2 
L\z 2 



+ Li) 



(d(*) 2 + C 2 (zf) 



,2 , V 2 * 2 



+ -^(Ciiz) 2 - C (z)C 2 (z)) + — - (p 2 z 2 + 3)d(^) 2 + ?—(C (z) + 3C 2 (z)) 



+ 



+ 



+ 



+ 



1 



-A^ + Li) 



L 2 (L 2 -L 2 ) 2 (AD-BC) 2 
L\z 2 



{C x {z? + C 2 {z) 2 ) 



-(C 1 (^) 2 -Co(^)C 2 (^)) + 



10p 2 



(pV + i)C x {z) 2 + ^(C (z) + 3C 2 (z))' 



-C 2 (z)C 2 (z) 



L 2 {L\ - L 2 ){AD - BC){AD - BC) \ 6 
j ~ ^f) C x {z)C x {z) + ^f(C (z)C 2 (z) + C 2 (z)C (z)) 
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Writing Ey(p) =Y 1 + Y 2 + Y 3 + Y 4 - 2Y 5 - 2Y 6 , we obtain 



Y 1 = 
Y 2 = 
+ 



3%M^L 



4L 2 {AD - BC) 2 
3iM 5 Lp 2 



z 2 CAz) 



Li 
Lo 



— lIc (l ) 2 , 



T 2 - T 2 



4L 2 (L 2 - L 2 )%AD - BC) 2 \ p 8 
p 6 z 4 L 2 



p 6 z 6 



[C^C^ + C^C^z)] 



[C^zf + C 2 (z) 2 ] - L\p'z 2 \C x {z) 2 - C (z)C 2 (z)} 



4 4 

p Z 



2p 2 L{ [ CoWtf^aiz) - 4p 2 z 2 C 2 (z)] + ^[ Co (z)C 2 (z) + Ci(*)C 3 (*)] 



Li 



+ - [(L 2 F(L ) - L?F(L 1 ))(G(L 1 ) - G(L )) - (LgG(L0 - L?G(L ))(F(L 1 ) - F(L ))} 



+ " ^o)) 2 + ^(F(Li) ~ F(L )) 2 ) , 



Y, 



3iM 5 L 



p A Li{Li-L*)C Q {Li)C*{L x ) 



4L 2 (L 2 -L 2 )(AD-BC) 2 p* 
+ p 4 Ll{L\ - L 2 )[C (L ) 2 + C^Lo) 2 } + 4p 2 ^(A,) 2 ' 



+ 



3iM 5 Lp 2 



AL 2 (L 2 - L 2 ) 2 (AD - BC) 2 
L 2 z 4 



■ — [C 1 (z)C 3 (z)+C 2 (z)C^z)] 
op 



[Ci(z)(C 3 (z) - Ci(z)) + C 2 (z)(C (z) - C 2 (z))] 



- ^C'o(L 1 )C' 2 (L 1 ) 



3p 2 

2^0 ) (L 1 )C 3 (L 1 ) - ^[C^Lo) 2 + O,(L ) 2 ] - 2^f tf 2 (Z„)[Zr,< '*< L„) - - LjC ( U)] ) 



7-47-2 
-^1-^0 I 



3iM 5 Lp 2 



4L 2 {L 2 - L 2 ) 2 (AD - BC){AD - BC) 
z 2 L\ 



2 (L °~ L?) ^[^(^^(z) + (7 2 (z)C 2 (z)] 



+ 



4 [2^(^(7!^) -Ca^Co^)- ^0(^(72^)^ - 2 L ° C ^ Lo) [G(z) - p 2 L\F{z) 

AL 2 {L 2 -L 2 ){AD-BC){AD-BC) \ ? LAAW + L.C^L^ 

- & \^-^-C 2 (Li)C 2 (Li) + 2p 3 L?C 1 (L 1 )C (L 1 ) + J(8 +p 2 L 2 )p 2 L 2 C 2 (L 1 )C (L 1 ) 

1 



Li 



- |(8 + p 2 L 2 )p 2 L 2 [2C 1 (Lo)Ci(Lo) - C (L )(7 2 (Lo) - C 2 (L )C , (L )] 



„4 r4 



[d^Ci^o) + C 2 (L )C 2 (L )] - 2p 3 L3[C 1 (L )C' (L ) + C 2 (L )C 1 (L )] 
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C The axial two-point function 



The axial two-point function is defined by 

< j%{x)r A {x') >= -i J dV'^'^n^Go). 

The Ward identities require a transversal 11^ (p), so that we may write H A v (p) = Il^(p 2 )T Mly , 
where T^ u is the transversal projector in p. To calculate this two-point function in our 
AdS/QCD setup, we have to use the boundary Lagrangian as shown in Eq. (fT2]) . This allows 
us to write an expression for the axial current by differentiating the full Lagrangian with 
respect to the boundary source a M . Schematically, we obtain the following expression 

r A (x) = M*±[-d z Ap{x, z) + d,A 5 (x, z)\ . 

Z Lo 

Plugging this expression into the equation for the axial two-point function, we find that we 
can write the result as the sum of the propagators of the and A 5 fields, giving 

< J%{x)r A {x') > = ^d z ,(^d z <A tl {x,z)A u {x'z')>\ 

Z 1 V Z J z,z'=L 

dld v <A 5 (x,L )A 5 (x',L )>. 

Now, making use of the definitions of the axial propagators found in section 5 above, it is 
easy to see that the A 5 propagator cancels the longitudinal part of the A^ bulk-to-boundary 
propagator. Thus, the boundary term containing A 5 is essential for the satisfaction of the 
transversality condition. 
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D Simplification to transverse boundary propagators 



The proof that only the transverse part of the boundary propagators is necessary for the cal- 
culation of Wllrr{Q 2 ) and Wlrlr(Q 2 ) is as follows: Consider a general four-point function 
with external momenta p M , p u , l a and lp. We assume the most general boundary propagators, 
i.e. with transverse and longitudinal parts. We can write this general four-point function as 

+ d 1 T a T p T ll p v + d 2 T a T (3 p^T u + d z T a l p T^T u + d 4 l a T /3 T^T u 
+ e l T a T f5 p^p y + e 2 l a l (3 T^T u + e 3 T a l^p u + e A T a fp»T u 
+ eslTyT" + e & l a T?T»p u 

+ hT a fp»p v + f 2 l a T p p^p v + fol a fT»p v + f A l a fp^T v 
+ gi rfp»p u . 

In this expression, T M and T v are transverse projectors in p M and p v respectively, whereas T a 
and T 13 are transverse projectors in l a and I 13 . We omit the second index of every transverse 
projector, since it is contracted inside the coefficient functions q — > gi, where all the vertex 
and gauge structure resides. For instance, d l T a T (i T^p v is really d^p^T^'T^'T^ 'p v p v \ 
where this term can originate from the longitudinal part of a bulk-to-boundary vector or 
axial-vector propagator, or from an particle emitted from an axial-vector current at the 
boundary. Now, the Ward identities obeyed by the four-point functions calculated above are 

IJpWr^O and p,p„Wr^ = ■ 

Applying the first Ward identity to the general form for the correlator, we obtain that 

e 2 T^T u + 9l p»p v + ftffT v + / 3 TV = 0. 

But this can only mean that each term in this equation separately vanishes, because these 
terms are linearly independent. Arguing similarly, one can apply the second Ward identity 
to find that the coefficients which must be zero are ei, e 2 , j\, f 2 , f'3, and g\. Now, what 
we are really after is the Lorentz singlet quantity given by the contraction of W^ va ^ with 
ViivVap (the equivalent of Wllrr(Q 2 ) and Wlrlr{Q 2 ))- This contraction trivially removes all 
the remaining non-zero terms, apart from the C\ term, which is precisely the one composed 
entirely of transverse external propagators. Thus, we have shown that the only term that 
contributes to the scalar functions Wllrr{Q 2 ) and Wlrlr{Q 2 ) is the one obtained by using 
the transverse propagators only on the external lines of the Witten diagrams. 
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